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experimental data is performed. 
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1. Introduction 

The physics of Ught neutron-rich nuclei is particularly interesting, because of the appearance 
of their exotic structures, like those of halo nuclei. The lightest nuclei of this kind are found in 
the Helium isotope chain: ^He as a two-neutron and ^He as a four-neutron halo nucleus. They are 
both radioactive and undergo jS-decay with a half life of = 0.8s and 0. Is, respectively. Despite 
the short life time, a combination of atomic and nuclear physics techniques, has enabled precise 
measurements of ground-state observables like the energy and the charge radius [jl]]. Excited-state 
properties, like electromagnetic transitions in the continuum have been investigated in the past 
with Coulomb dissociation experiments by Aumann et al. [^. Tackling the theoretical study of 
these nuclei is very challenging, because one needs to simultaneously describe the small separation 
energy of the halo neutrons and the large radius of the whole system. Because they are light-mass 
nuclei, one can use ab-initio techniques to study them. Here, we will discuss the recent results 
obtained using the hyperspherical harmonics (HH) method. 

A major breakthrough in nuclear physics has been the development of chiral effective field 
theory, which is well routed to Quantum Chromo Dynamics. Even though several light nuclei have 
been investigated with chiral potentials, we are still missing a prediction of Helium halo nuclei from 
chiral Hamiltonians. Here, we show a first step taken in this direction by using low-momentum 
chiral two-nucleon forces. 

The paper is organized as follows. In section 2 we will introduce the hyperspherical harmonics 
method. In section 3 and 4 we will present results for ground-state and excited-state properties, 
respectively. Finally, in section 5 we will draw some conclusions. 

2. Hyperspherical Harmonics 

Given the Hamiltonian H we use the HH expansion to solve the Schrodinger equation. The 
HH method is typically a few-body method used for 3 and 4-body systems. Using the powerful 
antisymmetrization algorithm introduced in [||], it is possible to extend the method to a larger mass 
number and tackle HeUum halo nuclei. The HH approach starts from the Jacobi coordinates 

1 ^ ( 1 \ 

where r, are the particle coordinates. Using the t], one can then transform to hyperspherical coor- 
dinates composed of one hyperradial coordinate p = \Jy^^Ii J7? and a set of (3A — 4) angles that 
we denote with Q. (for more details see [^). Using this coordinates one can recursively construct 
the hyperspherical harmonics ^j^] and use them as a complete basis to expand the wave function. 
Such expansion reads 

*I'(r]i,...,T]^_l,5l,...,5A,fl,"-,fA) = XI L C[/f]„/?„(p)%](^l,5i,...,5A,fl,"-,fA), (2.2) 

n [K] 

where Si and tj we, the spin and isospin of the nucleon i, respectively; C[^]„ is the coefficient of 
the expansion, labeled by [K], which represents a cumulative quantum number that includes the 
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grandangular momentum K; n labels the hyperradial wave function R „{p). The latter is expanded 
in terms of the generalized Laguerre polynomials times an exponentially falling off function, which 
is essential to speed up the convergence of halo nuclei, due to their extended tail. To further 
increase the convergence rate of the calculations, we typically employ an effective interaction in 
the hyperspherical harmonics (EIHH), as first introduced in [Q]. 

3. Results for ground state properties 

In the following, we present our results for the ground state energy and the point-proton radius 
calculated using the hyperspherical harmonics method. As input Hamiltonian we employ a class of 
low-momentum potentials, which are obtained applying a Viow/t procedure [||] on a starting chiral 
two-body potential ^ at next-to-next-to-next-to leading order (N ^LO). 

In Fig. |l], we show the ^He convergence patterns of the EIHH method for the ground-state 
energy and the point-proton radius as a function of the maximal grandangular momentum [Q]. 
Three different cutoffs A = 1.8,2.0,2.4 fm^^ of the Viowk chiral potential are shown. One observes 
that the convergence is very nice for all A's, even for the largest cutoff. Concerning the energy we 
also show that the extrapolated results from previous work on hyperspherical harmonics [Q| agree 
nicely with EIHH. The convergence rate is very good also for the point-proton radius, which allows 
us to provide solid results for ^He. 

We have also explored the four-neutron halo nucleus ^He with hyperspherical harmonics. In 
Fig. ^, we show the ^He ground state energy from a Viowk chiral potential with A = 1.8 fm^' as 
a function of K^ax- We present both the variational HH expansion (where we do not apply the 
effective interaction) and the EIHH results. The convergence for ^He is quite slow, as indicated 
by the fact that the HH and EIHH patterns do not merge yet at ^max = 10. Interestingly, an ex- 
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Figure 2: (Color online) The He ground-state energy calculated with the HH and EIHH methods as a 
function of obtained for a cutoff of A = 1.8 fm^' of the Viowjt chiral potential. The extrapolated HH 
results are shown as a reference and compared to the Coupled-Cluster results. 



2.1 



IS 



1.7 



' ' I 

• GFMC 
□ NCSM 
vFMD 
AMCM 

• EIHH 



0.2 




0.4 



0.6 



0.8 



[MeV] 



■■■■*"^lowk(N"LO) 



1.2 



1.4 



Figure 3: (Color online) Correlation plot of the ^He point-proton radius versus the two-neutron separation 
energy 5'2„. The experimental range is compared to theory based on different ab-initio methods (see text). 



trapolation of the variational HH data lies close to the Coupled-Cluster result from with the 
same interaction. Nevertheless, at the moment we are not able to provide precise results for ^He 
from hyperspherical harmonics. Thus, we concentrate on analyzing the correlation between en- 
ergy and radius just for ^He. In Fig. ^ we plot versus the two-neutron separation energy S 2n 
and present a combined comparison of our results to experiment and other ab-initio calculations: 
Green's Function Monte Carlo (GFMC), No Core Shell Model (NCSM), Fermionic Molecular Dy- 
namics (FMD) and Microscopic Cluster Model (MCM) (see also |^ and references therein). The 
cutoff dependence of our results with Viowi allows us to study the correlation between these observ- 
ables: the radius increases as the separation energy decreases. Our calculations do not reproduce 
simultaneously r p and S2n '■ there exists an optimal value of A where S2n is predicted in accordance 
with experiment, but is not reproduced and vice-versa. Also, we would like to note that other 
calculations which omit JiNF, (all except from the GFMC) do not go though the experimental band. 
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Figure 4: (Color online) Panel (a): The coiTelation between «£ and Sin in ^He obtained with the the 
Minnesota potential varying the parameter u. Panel (b): The correlation between rs^n and Ue in different 
model spaces (different Xmax)- Shown is also the extrapolated value of rskin, which is used to estimate Ue- 



This points towards the importance of including three-nucleon forces in the Hamiltonian. 

4. Results for excited state properties 

As an example of excited-state properties of halo nuclei we report about our recent calculation 
of the electric dipole polarizability of ^He [Q]. is related to the inelastic response of the 
nucleus to an externally applied electric field and is relevant in the extraction of nuclear quantities 
from atomic spectroscopic measurements. The atomic energy levels, in fact, are affected by po- 
larization of the nucleus due to the electric field of the surrounding electrons. The polarizability 
of ^He could be extracted from Coulomb dissociation measurement of the dipole transition by 
Aumann et al. ^ and was reported in Ref. [ |To| ] to be much bigger than the polarizability of '^He. 
The electric dipole polarizability is defined by 

where I^I'o//) is the ground state and final state of the nucleus and £1 is the dipole operator. Because 
it requires the knowledge of the dipole spectrum of the nucleus, its theoretical evaluation is more 
involved than a bound-state calculation. We perform our calculation with the EIHH method by 
using the Lanczos algorithm with a starting dipole pivot as explained in [^. As nuclear potential 
we chose the simple semi-realistic Minnesota force, which reproduces the experimental value of 
the polarizability of '^He reasonably well. Within this force model we can add attractive P— wave 
interactions by changing the parameter u (see for details). This mostly affects ^He, without 
substantially changing "^He. By varying u we first observe a correlation of Ue vs S2n, as shown in 
Fig. ^(a). We have chosen u so that the halo feature, represented by S2n, is reproduced. We then 
study the correlations between and the skin radius rgkin = rn — rp, where r„ is the mean point- 
neutron radius. By varying the model space we observed that and rgidn are correlated linearly 
for ^max > 6 as «£ = a + rgkin- From our theoretical data we fit the coefficients a and b and then 
we used them to estimate the polarizability out of a bound-state calculation of the skin radius. The 
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calculation of r skin, in fact, does not require an expansion on the dipole excited states and as such 
is less computationally demanding and can be performed for larger model spaces {K max = 16) and 
then extrapolated exponentially, leading to rgkin = 0.87(5) fm. Using our extrapolated skin radius 
and the linear dependence, we estimate the theoretical nuclear electric polarizability of ^He to 
be = 1.00(14) fm^. The error bar is obtained by propagating the errors on a, b and rskin- We 
observe that our theoretical estimate is about a factor of two smaller than the experimental value 
of a^"^ = 1.99(40) fm^ [10|. This points toward a potential disagreement between theory and 
experiment. Investigations with chiral potentials can possibly help understanding this discrepancy. 



5. Conclusions 



In conclusion, we have presented our recent results on Helium halo nuclei from hyperspheri- 
cal harmonics. The binding energy and the radius can be precisely calculated for ^He using chiral 
low-momentum two-body forces. The obtained cutoff dependence together with a comparison to 
the experiment serves to highlight the importance of three-nucleon forces. We also discussed the 
nuclear dipole polarizability of ^He as an excited state observable, which has recently attracted at- 
tention. Our estimate from simplified nuclear potentials leads to a disagreement with experimental 
data, which will be hopefully clarified in the future when realistic chiral potentials will be used. 
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